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Abstract 

"Causal Dynamical Triangulations" (CDT) represent a lattice regularization of the 
sum over spacetime histories, providing us with a non-perturbative formulation of 
quantum gravity. The ultraviolet fixed points of the lattice theory can be used to 
define a continuum quantum field theory, potentially making contact with quantum 
gravity defined via asymptotic safety. We describe the formalism of CDT, its phase 
diagram, and the quantum geometries emerging from it. We also argue that the 
formalism should be able to describe a more general class of quantum-gravitational 
models of Horava-Lifshitz type. 
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1 Introduction 



At this stage, there is no certainty how to best reconcile the classical theory of rel- 
ativity with quantum mechanics. Applying the well-tested methods of quantization 
to gravity - defined by the Einstein-Hilbert action - and quantizing the fluctuations 
around a classical solution to Einstein's equations leads to a non-renormalizable 
theory. This happens because in four spacetime dimensions the mass dimension of 
the gravitational coupling constant G (in units where h and c are 1) is —2, whereas 
it should be larger than or equal to for the theory to be renormalizable perturba- 
tively. One would therefore expect the perturbative effective quantum field theory 
description to break down at energies E satisfying GE'^ > 1. 

There are of course well-known examples where the non-renormalizability of a 
quantum field theory in the ultraviolet (UV) was eventually resolved by introducing 
new degrees of freedom, missed initially because they were not directly observable 
at low energies. The electroweak theory is an example where perturbative renor- 
malizability was "regained" in this way. The theory was described first by a four- 
fermion interaction with an associated Fermi coupling Gp of mass dimension —2, 
just like the Newton constant G in gravity. As a result, its perturbation theory 
breaks down at energies with GpE^ > 1. However, it turns out that for energies 
above l/y/Gp ~ Mw, the mass of the VF-particle, the four-fermion theory has to be 
replaced by the S'C/(2)-gauge theory of the weak interactions, which contains new 
excitations, the W- and Z-bosons. The new electroweak theory is a renormalizable 
quantum field theory. 

Similarly, in the 1960s the low-energy scattering of pions was described by a 
non-linear sigma model, another non-renormalizable quantum field theory whose 
coupling constant, the pion decay constant F^r-squared, has mass dimension —2. 
However, high-energy scattering at energies beyond 1/-F^ is no longer described well 
by the non-linear sigma model, because it starts probing the intrinsic structure 
of the pions. A correct description has to incorporate appropriate new degrees of 
freedom, the quarks and gluons, and the corresponding quantum theory - quantum 
chromodynamics - is perfectly renormalizable. 

There is no obvious reason which prevents us from writing down a perturbative 
(and non-renormalizable) expansion for gravity around some classical background 
geometry, say, fiat Minkowski spacetime, if we are interested in an effective quan- 
tum field-theoretic description whose range of applicability does not extend beyond 
energies with GE'^ fa 1. In view of the examples cited above, it is then tempting to 
conjecture that the apparent non-renormalizability of gravity could be resolved by 
the appearance of new degrees of freedom at higher energies, rendering the theory 
renormalizable after all. 

A solution of this kind may be in the form of a superstring theory in a higher- 
dimensional spacetime, where the gravitational excitations are intertwined with in- 
finitely many new degrees of freedom in such a way as to cure the UV problem. 
Although string theory cannot be ruled out as the correct answer, the world picture 
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it provides has yet to be verified. In particular, supersymmetry - predicted by string 
theory - has not yet been observed at the Large Hadron CoUider. Of course, even 
if no evidence of supersymmetry is found at this or future coUiders, it may still be 
present at even higher energies. In this sense, the absence of observational evidence 
for supersymmetry does not disprove superstring theory as such, although it makes 
it less compelling as a resolution of the problem of unifying gravity and quantum 
theory. 

There are other potential resolutions to the problem of finding a suitable "ultravi- 
olet completion" of perturbative quantum gravity, which are not based on fundamen- 
tal, string-like excitations and do not obviously require the existence of supersym- 
metry or extra dimensions. These are so-called non-perturbative approaches, whose 
starting point typically consists of a set of dynamical degrees of freedom closely 
modeled on those of classical gravity ("curved geometry" in one way or other), to- 
gether with a non-perturbative prescription for quantization. A concrete example, 
that of Causal Dynamical Triangulations, will be described in some detail below. 
Its geometric degrees of freedom, in presence of a UV cut-off, are given in terms 
of triangulated, piecewise flat spacetimes with discrete curvature assignments. Its 
non-perturbative quantization follows that of a standard lattice field theory, albeit 
with a dynamical rather than a fixed lattice. 

An obvious charm of such a purely quantum field-theoretic ansatz lies in its 
minimalism, and the absence - to a large degree - of free parameters and other 
"tunable" ingredients. On the other hand, a key difficulty of this type of approach 
is to demonstrate that it is related to classical gravity in a suitable limit, something 
that is not at all obvious once one has moved beyond linearized quantum fields on 
a fixed background spacetime. One also needs to spell out what it means for the 
non-perturbative theory to exist, which likewise is non-trivial in a background-free 
description where "observables" are hard to come by. 

In parallel with advances in string theory, also research in the wider area of non- 
perturbative quantum gravity has seen a steady rise in interest in recent decades. On 
the one hand, this was due to the rejuvenation of canonical quantum gravity in the 
form of loop quantum gravity from the late 1980s onwards.^ At about the same time, 
the covariant gravitational path integral was given a new, non-perturbative lease of 
life in terms of "dynamical triangulations" . Motivated originally by the search for 
a non-perturbative dynamics of curved, two-dimensional worldsheets in (bosonic) 
string theory, this dynamical lattice formulation provides a powerful computational 
tool for evaluating gravitational path integrals quantitatively: analytically in two, 
and numerically in higher dimensions. - The focus of the present article will be 
on this latter development, arguably the conceptually most straightforward and 
methodologically minimalist extension of the standard perturbative and covariant 
quantum field-theoretic formulation of gravity. We will explain how it may lead to 
the construction of a viable theory of quantum gravity, valid on all scales, without 

^Curiously, this ansatz also postulates the fundamental character of certain one-dimensional 
"closed-string" (a.k.a. "loop") excitations in the quantum theory. 
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running into contradictions vis-a-vis the perturbative non-renormalizability of the 
theory. 

In the late 1970s, Weinberg outhned a scenario, coined asymptotic safety [1], for 
how quantum field theories which are not power-counting renormalizable around a 
trivial Gaussian fixed point could under certain, general conditions still make sense, 
just like ordinary renormalizable theories. In particular, an asymptotically safe 
theory is characterized by only a finite number of coupling constants, whose values 
will be determined by comparison with experiment or observation. The asymptotic 
freedom scenario is naturally described in the language of quantum field theory 
and the renormalization group. It is characterized by the presence of an ultraviolet 
fixed point in the infinite-dimensional coupling constant space of a theory, with the 
property that in the fixed point's neighbourhood the dimension of the subspace 
of attraction is infinite-dimensional, with finite co-dimension. This co-dimension 
coincides with the number of free parameters of the theory that need to be fixed by 
experiment. Such a UV fixed point therefore attains a similar status to that of the 
Gaussian fixed point of a renormalizable theory. The snag is that the tools of the 
perturbative theory are usually not sufficient to find such ultraviolet fixed points - 
if they exist for a given theory - and to study their neighbourhoods. 

To illustrate the implications of the presence of such a fixed point (in a somewhat 
simplistic fashion), let us introduce the dimensionless coupling 

G{E) := GE^. (1) 

A fixed point in this context always refers to the behaviour under a change of scale E 
of a dimensionless, energy- dependent function like G{E). The dimensionful quantity 
G in (1) can at this stage still be thought of as a (classical, low-energy) coupling 
constant of mass dimension —2. Let the behaviour of G{E) be dictated by a beta 
function (5{G) according to 

E— = ^{G), with ^{G) = 2(7 - 2u&, (2) 
(xE 

for some real parameter u. It is immediately clear that for a; ^ 0, G = const 
is no longer a solution to (2). For consistency, G has to acquire a non-trivial E- 
dependence and therefore becomes a function G{E) — G{E)/E'^. In (2) we have 
chosen the simplest non-trivial beta function such that (i) in the limit of low energy, 
E ^ 0, G{E) goes to a constant (which we will continue to call G), and (ii) for 
E ^ oo, G{E) goes to a non-trivial UV fixed point. Explicitly, the solution to the 
differential equation in (2) can be stated as 

from which we can read off the location of the UV fixed point at G = the 
non-trivial zero of the beta function. An important feature of this solution is that 
the couphng constant G{E) goes to zero at the UV fixed point. 
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Figure 1: Changing an asymptotically free theory to an asymptotically safe one by 
increasing its dimension from dio d-\-e results in a shift of its ultraviolet fixed point 
to a value 5^ > 0. 

In case the above example should appear somewhat ad hoc, it can be understood 
as arising from a more general construction, which starts from an asymptotically 
free theory in d dimensions. Fig. 1 (left) illustrates the corresponding (negative) 
beta function of the coupling together with a Gaussian UV fixed point at (7 = 0. 
If this theory is "lifted" to d + e dimensions - assuming that such a perturbation in 
the dimension is well defined, at least for small £ > - its beta function will change 
according to 

P{g)^p{e)g + P{g), (4) 

where p{e) is the (positive) amount by which the mass dimension of the coupling 
g decreases as a result of the dimensional increase by e. (Our previous example, 
whose beta function was defined in relation (2), corresponds to p = 2.) Note that 
the Gaussian UV fixed point of the original theory has become a non-trivial UV 
fixed point away from zero in the higher- dimensional theory, while g = Q has been 
turned into an infrared fixed point, as illustrated by Fig. 1. 

The theories we have discussed so far - four-Fermi theory, non-hnear sigma 
model and Einstein gravity - display a similar behaviour in the sense that they are 
asymptotically free, renormalizablc theories in spacetimc dimension d = 2. Trying to 
make sense of them beyond dimension 2 by way of a 2-|-£-expansion, one encounters 
the situation depicted in Fig. 1. Of course, one may formally set £ = 2 in such an 
expansion, as would be needed to reach the dimension d — A oi physical spacetime, 
but the validity of the pcrturbative expansion for such large values of £ would need 
to be established to take the results seriously, and a priori appears perhaps rather 
doubtful. 

Non-trivial UV-complete extensions to 0? = 4 of the four-Fermi interaction or the 
non- linear sigma model are not known and presumably do not exist. As mentioned 
above, we should rather think of them as effective theories, which happen to describe 
certain low-energy properties of more fundamental theories with more and different 
fundamental excitations. Still, it is difficult to draw any conclusions from this for 
general relativity, the theory we are interested in, which is after all very different 
physically: exactly the degrees of freedom that are fixed in all other theories, those 
of spacetime itself, become dynamical in gravity. Much work has gone into trying 
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to show that four-dimensional gravity possesses an ultraviolet fixed point with the 
requisite properties, either in terms of the 2 + £-expansion [2] or by using general 
renormalization group techniques [3]. 

In what follows, we will not be concerned with the details of these efforts, but 
with the question of how the hypothesis of asymptotically safe gravity may be tested 
independently and non-perturbatively by using standard field-theoretic tools and by 
formulating quantum gravity via a lattice regularization. 

2 A lattice theory for gravity 

A number of issues have to be addressed when representing gravity on a lattice. Is 
it possible in principle to construct a well-defined lattice regularization of gravity 
with a UV lattice cut-off, which can be removed in a controlled way to obtain a 
continuum limit (whatever this may turn out to be)? The answer is yes. More 
precisely, the issue is not so much how to represent gravity on a lattice, but how 
to represent a theory as a lattice theory whose standard continuum formulation 
in terms of local fields is diffeomorphism- invariant, a vast gauge invariance closely 
related to the differentiable structure of the underlying manifold and its description 
in terms of local coordinate charts. 

For the geometric degrees of freedom of the gravitational theory this can be done 
by viewing the lattice itself as representing directly a (piecewise linear) geometry. 
The key point is that such a geometry can be described uniquely without ever in- 
troducing coordinates, thus circumventing the associated redundancy of having to 
choose any particular set of coordinates. A convenient choice is to use lattices which 
are triangulations, in the sense of consisting of rf-simpliccs, triangular building blocks 
which are ci- dimensional generalizations of flat triangles (=2-simplices). Assuming 
the interior of a (i-simplex to be flat, its geometry is uniquely specifled by giving the 
lengths of its d{d + l)/2 one-dimensional edges or links. Together with the infor- 
mation of how the simphces are "glued together" (that is, how {d — 1) -dimensional 
boundary simplices are identifled pairwise) to form a triangulated manifold, this suf- 
flces to compute all geometric information, including distances, geodesies, volumes 
etc. without using coordinates. Important for our path integral representation, 
Regge observed that the curvature of such a piecewise linear geometry is in a natu- 
ral way located on its (c^ — 2)-dimensional subsimphces (the "hinges"). By the same 
token, the scalar curvature term of the Einstein action of such a geometry is given 
by the sum over all hinges of the deficit angle around each hinge, multiplied by the 
hinge's volume [4]. 

In our construction of a theory of quantum gravity, the lattice-regularized path 
integral over geometries thus becomes the sum over such triangulations, with weight 
depending on the Regge implementation of the Einstein action. Precisely which 
class of triangulations should we sum over in the path integral? When applying 
Regge calculus to classical gravity one uses a fixed lattice, in the sense of leaving 
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the connectivity of its constituent simplicial building blocks unaltered. This still 
allows the curvature of the triangulation to be changed - for example, to optimally 
approximate that of a given smooth geometry - by changing the lengths of its one- 
dimensional edges. 

When using the piecewise linear geometries in a path integral, the task is dif- 
ferent. Firstly, we do not expect the individual path integral configurations to be 
smooth, but only continuous, in the same way as the paths in the path integral of a 
quantum-mechanical particle arc continuous but in general non-smooth^. Similarly, 
the piecewise linear geometries are a subset of all continuous spacetime geometries. 
Note that we can even restrict ourselves to a subset of piecewise linear geometries 
as long as it is suitably dense in the set of all geometries. More precisely, when the 
lattice spacing goes to zero, we require the expectation values of observables, again 
suitably defined on the piecewise linear geometries, to converge to the value they 
would take in the continuum quantum field theory (which we assume exists). In 
contrast with the aim of the classical theory, we are therefore not trying to approx- 
imate any particular geometry by our lattice geometries, but to span the whole set 
of geometries. 

In this context a specific subset of piecewise linear geometries has proved to be 
very useful, namely, the triangulations whose edges have all the same length, a, 
say. One can characterize this set of geometries as being constructed from gluing 
together equilateral simplicial building blocks in all possible ways, compatible with 
certain constraints (typically, a fixed topology and fixed boundary components). 
Consequently, the variation in geometry (the way in which the geometric degrees 
of freedom are encoded) is linked to the mutual connectivity of the building blocks 
created by the gluing and not to variations in the link lengths, giving rise to the 
name Dynamical Triangulations (DT) [5, 6, 7]. From a path-integral perspective this 
approach has the advantage that distinct triangulations correspond to physically dis- 
tinct geometries. Summing over this DT ensemble of geometries may therefore lead 
directly to the correct continuum measure in the limit that the UV cut-off is taken 
to zero, a ^ 0. By contrast, treating the triangulations classically a la Regge, with 
fixed lattice connectivity and variable link lengths, still contains redundancies, in the 
sense that many different lattice configurations can correspond to the same physical 
geometry (see [8] and references therein). For illustration, consider a rectangle in 
the two-dimensional plane and triangulate its interior. Clearly, the interior vertices 
can be moved around locally in the plane without changing the flat geometry of 
the rectangle. However, since all of these are different as Regge triangulations, this 
leads to a severe overcounting in the path integral of quantum Regge calculus, for 
which there is currently no known fix. 

Most importantly, the viability of the DT lattice regularization has already been 
demonstrated in a non-trivial case, that of gravity (coupled to matter) in two dimen- 
sions. As mentioned above, two-dimensional gravity is a renormalizable quantum 

^in fact, with unit probability they are nowhere differentiable 
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field theory and various observables can be calculated analytically [9] . The dynam- 
ically triangulated two-dimensional lattice theory can also be solved, a number of 

observables can be calculated analytically and its continuum limit, taking the lattice 
spacing a — ?> 0, can be taken [10]. Remarkably, results from the two different calcu- 
lations can be compared and are found to agree. We conclude that it is possible to 
provide a viable lattice regularization of a dijfeomorphism-invariant quantum theory 
of geometries. 

One may object that this two-dimensional theory has little to do with true grav- 
ity in four spacetime dimensions; to start with, it has no propagating gravitons. 
However, we would like to argue that it is much more a theory of fluctuating geome- 
tries than one would ever expect of the four-dimensional theory. Because there is 
no Einstein- Hilbert action in two dimensions (it is topological), each configuration 
contributes in the path integral with the same weight, which is a maximally quan- 
tum situation. This is borne out by the analytic solutions of this model, which show 
the two-dimensional geometries as wildly quantum-fluctuating. Nevertheless the 
lattice theory has no problem in reproducing the correct diffeomorphism-invariant 
continuum theory, also known as quantum Liouville gravity. 

2.1 Observables 

How to define what does and does not constitute an "observable" in quantum gravity, 
and how to construct and evaluate observables in any given formulation are physical 
questions of central importance . What we would like to highlight here is that a 
beautiful aspect of a geometric lattice formulation of quantum gravity of the type 
we are considering is that it forces one to address such questions head-on. It is not 
possible to hide behind some "expansion around fiat spacetime" , but one is forced 
to think in terms of physical "rods and clocks", much in the spirit of Einstein's 
classical theory. 

Let us discuss the basic objects of any quantum field theory, namely, the correla- 
tors of local quantum operators 0{x). Such correlators are important ingredients in 
constructing S-matrix elements, i.e. observables in quantum field theory on a fixed 
background. Also in conventional lattice theories, correlators play a crucial role in 
showing that a lattice theory has a continuum limit when the lattice spacing goes 
to zero. 

Consider some lattice scalar field theory, and let 0{xn) be an operator at lat- 
tice spacetime coordinate Xn — n ■ a, where a is the lattice spacing and n the 
integer-valued lattice coordinate. In general, we expect the correlator to fall off 
exponentially, 

- log(C(x„)C(x^)) ~ |n - m\/^igo) + o(|n - m|), (5) 

where go is the bare lattice coupling and C,{go) the correlation length in lattice 
spacings. The standard procedure for a lattice system is to take the continuum limit 
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at a second-order phase transition point g^, where the correlation length diverges 
hke ^ 

^(^o) oc — , a(^o) oc 1^0 - 9o\''- (6) 

Ifi'o — fi'ol 

Eq. (6) tells us at what rate we should scale the lattice spacing to zero in the limit 
go — >■ gQ, in order to find an exponential decay in the continuum, when the lattice 
correlation diverges, but the (dimensional) physical length Xn — Xm — {n — m)a is 
kept constant, 

mpuaigo) = l/i{9o), e-l"-"^l/«(^°) = ^-m,u\xr.-x^\ _ 

Eq. (7) illustrates the fact that dimensionful observables, like the physical mass niph, 
are defined by the approach to the critical point, not at the critical point. 

The existence of a critical point and an associated divergent correlation length 
constitute the backbone of the Wilsonian renormalization group approach to quan- 
tum field theory. Since we are appealing to this Wilsonian approach by asking 
whether asymptotic safety is realized, it is important to understand whether it 
can be applied to quantum gravity at all. A first step in this direction is to un- 
derstand whether suitable correlators and a correlation length can be defined in a 
diffeomorphism-invariant theory like quantum gravity. To start with, how can we 
define the distance between two points in a path integral where we integrate over 
the geometries defining this distance? 

In fiat (i- dimensional spacetime, let us rewrite the correlator of a scalar field 
say, in the form 

{H{R))v = ^ ^ y P^e-^I-^] jd'x Jd'y 6{R-\x - y\). (8) 

As indicated, this expression depends on a chosen distance R, but no longer on 

specific points x and y, which instead are integrated over. The integrand can be 
read "from right to left" as first averaging over all points y at a distance R from 
some fixed point x, normalized by the volume s{R) of the spherical shell of radius R, 
and then averaging over all points x, normalized by the total volume V of spacetime. 
We assume translational and rotational invariance of the theory and that V is so 
large that we can ignore any boundary effects related to a finite volume. 

This definition of a correlator is of course non-local, but unlike the underlying 
locally defined correlator has a straightforward diffeomorphism-invariant generaliza- 
tion to the case where gravity is dynamical, namely, 

mmv = ^jns\ fvi,f4,e-''^-i'U(v-jd'xVd^)- 
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which now includes a functional integration over geometries^ [g], and dependences 
of the action, measures, distances and volumes on [g]. Can the definition (9) be 

implemented meaningfully to define correlators in a quantum gravity theory? The 
answer is yes, and a two-dimensional example can again be used to demonstrate 
this. Namely, there are analytic predictions for the behaviour of the propagators 
of certain matter theories coupled to two-dimensional Euclidean gravity [9], which 
have been shown to be reproduced by numerical simulations of the corresponding 
lattice theory [11]. By the way, their behaviour is quite different from that of the flat 
space correlators, another manifestation of the fact that two-dimensional gravity is 
a theory of strong geometric fluctuations. 

2.2 Time-slicing and baby universes 

An interesting aspect that can be analyzed in detail in the solvable two-dimensional 
quantum theory of fluctuating geometry is that of proper time . One usually con- 
siders a situation where the rotation to Euclidean signature has taken place and 
"proper time" is simply given by "geodesic distance". In this setting, a closed 
one-dimensional spatial universe of fixed "time" is simply a loop of length £. In the 
corresponding quantum theory one can ask for the amplitude for a universe of length 
£i to "propagate" to another one of length £2 in proper time t. More precisely, the 
outgoing loop of length £2 is said to have a proper-time (in this case a geodesic) 
distance t to the incoming loop of length £1 if each point on £2 has geodesic distance 
t to £1. (The geodesic distance from a point to a set of points is defined as the 
minimum of the geodesic distances from the point to the points in the set.) 

Fig. 2 shows a typical geometry in the path integral contributing to the corre- 
sponding amplitude G{£i, £2; t). It will often be convenient to work with its Laplace 
transform, 

/»00 /"OO 

G{x,y;t)^ / / d£id£2e-''^'-y^' G{£ij2;t). (10) 
Jo Jo 

We can view x and y in this expression as boundary cosmological constants, since x ■ £ 
would be the action of a one-dimensional "spacetime" of volume £ and cosmological 
constant x. 

As shown in [12], the amplitude G{x,y;t) satisfies the remarkably simple equa- 
tion 

dG{x,y,t) _ d{W{x)G{x,y,t)) 

dt dx ' ^ ' 

where W{x) is the Hartle-Hawking disk amplitude, which in two-dimensional Eu- 
clidean gravity is given by [10] 

W{x) = {x-'^)\lx + ^/lv. (12) 

^in accordance with standard notation, [g] denotes an equivalence class of metrics g under the 
action of the diffeomorphism group 
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Figure 2: Incoming and outgoing boundary loops of length £i and £2, separated by a 
geodesic distance t, and a typical interpolating geometry of cylinder topology which 
contributes to the amplitude £2; t) in Euclidean signature. The additional loops 
drawn onto the interior geometry consist of points which share the same distance to 
the incoming loop. As indicated by the upper set of three loops, there can be many 
disconnected loops at a given distance to the incoming loop. 

As is clear from Fig. 2, space can branch out into many disconnected parts (i.e. 
change its topology) as a function of proper time t, giving rise to "baby universes" . 
The appearance of baby universes on all scales leads to the two-dimensional quantum 
spacetime being fractal , with Hausdorff dimension dh = 4: [12, 13]. 

Rather amazingly, it is possible to integrate analytically over these baby uni- 
verses, resulting (for each time history) in a spacetime with a proper-time foliation 
and no baby universes [14]. Alternatively, the expression for the loop-loop prop- 
agator without baby universes can be obtained directly by summing over a class 
of two-dimensional spacetimes which from the outset lack baby universes, provided 
one redefines the coupling constants suitably [15]. This latter procedure can be 
implemented also at the regularized level in terms of a set of causal dynamical tri- 
angulations" (CDT), to be distinguished from the larger class of merely "dynamical 
triangulations" (DT), which served as carrier space for the Euclidean gravitational 
path integral [15]. 

The resulting theory has a well-defined Hamiltonian and corresponding unitary 
proper-time evolution. The explicit map between the cosmological constants of DT 
and CDT turns out to be non-analytic. 



where we have denoted the CDT-analogues of the couplings with a subscript and 
tilde. Consequently, in CDT both lengths and areas acquire a dimensionahty dif- 





(13) 
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Figure 3: The light cone structure (and therefore the underlying Lorentzian geom- 
etry) becomes degenerate in points where space splits in two. 

ferent from that found in the DT ensemble of spacetimes and in Liouville gravity. 
When using the CDT ensemble, also the Hausdorff dimension changes from 4 to 2, 
the canonical value for ordinary smooth two-dimensional spacetimes.^ 
The CDT loop-loop propagator satisfies the equation 

dG{x,y,t) _ d{{x^ - Acdt)Gix,y,t)) 

dt di ' ^ ^ 

and the Hamiltonian governing the (proper-) time evolution is given by 

Gih,i2, t) = (£2|e-*^|£i), H = -fi + K^tl (15) 

while the CDT Hartle-Hawking wave function (which is derived from the propagator 
G [15]) satisfies 

HW^dt{i) = 0. (Wheeler-DeWitt) (16) 

Above, our first way of deriving this formulation was as a kind of "effective" 
theory: we started from the set of all Euclidean two-dimensional geometries of a 
fixed topology. These geometries are "isotropic" in the sense that they do no carry 
any a priori preferred direction. We then superimposed a notion of proper time on 
them and integrated out part of the degrees of freedom. However, when starting 
in the physically correct Lorentzian signature, one can formulate a general principle 
which excludes geometries whose spatial topology is not constant in time [16]. The 
point is that spatial topology changes are associated with causality violations of one 
kind or other. This is illustrated by the "trouser geometry" depicted in Fig. 3. As is 
clear from the embedding of this two-dimensional spacetime in fiat Minkowski space, 
with time pointing upward, there must be at least one point near the crotch of the 
trousers where the tangent plane is exactly horizontal and the light cone therefore 

word of warning: the coincidence in Hausdorff dimension does not allow one to conclude 
that the quantum geometry of two-dimensional CDT in any way approximates a smooth classical 
manifold; in fact, it does not. 



12 



degenerate. Note that imposing causality conditions on the geometry to ehminate 
such configurations only makes sense in the presence of a Lorentzian metric and 
cannot even be formulated in a purely Euclidean theory, in the absence of any extra 
structure. 

By the same token, one can take as domain of the path integral the set of all 
Lorentzian piecewise fiat triangulations whose causal structure is well defined, and 
where in particular no changes of spatial topology are allowed to occur. The set of 
causal dynamical triangulations (CDT) - which can be defined in any dimension (not 
just d = 2) - obeys a strong version of causality of this kind, which is implemented 
by requiring each triangulation to be the product of a one-dimensional "triangula- 
tion" (a line with equidistant points), representing discrete proper time, and other 
triangulated degrees of freedom, representing the spatial directions of the geometry, 
which may be thought of as triangulated fibres over a one- dimensional base space. ^ 
As an added bonus, each triangulation in the class of CDT can be analytically 
continued to Euclidean signature, and the associated gravitational Regge actions 
satisfy the standard relation between actions defined in spacetimes of Lorentzian 
and Euclidean signatures, namely, 

Lorentzian ' ^ S j^^f^Hfi^an' (■^'^) 

Despite the fact that the actions obey (17), the Lorentzian theory defined on CDT 
geometries will even after this "Wick rotation" be distinct from the full Euclidean 
theory, because not every Euclidean triangulation is the image of a causal, Lorentzian 
one. The subclass of Euchdean geometries that are in the image can be obtained 
"surgically" as explained above, by superimposing a notion of proper time on each 
Euclidean triangulation and then removing all of its baby universes associated with 
spatial topology changes. The two-dimensional case is sufficiently simple to allow 
us to perform the calculation in either way, by starting from a path integral over 
all Euclidean geometries and removing baby universes, or by starting from a path 
integral over causal (CDT) geometries and rotating it to Euclidean signature. Both 
results agree after a redefinition of the coupling constants. Let us note in passing 
that our formulation - not only in dimension 2, but also in higher dimensions - has 
a couple of characteristics reminiscent of so-called Hof ava-Lifshitz gravity , namely, 
the use of a preferred time foliation and a unitary time evolution. We will return to 
this subject in Sec. 4 below. 

2.3 CDT in higher dimensions 

It is not known whether the above-described procedure of integrating out baby uni- 
verses in d = 2 can be generalized to higher dimensions in a simple and useful way. 
It implies that at this stage we have two a priori unrelated lattice gravity theo- 
ries in dimension d > 2, one purely Euclidean based on DT and one Lorentzian 

^Product triangulations, of which this is a particular instance, were investigated in [17], see also 
[18]. 
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based on CDT. The latter starts out in physical, Lorentzian signature, and imposes 
local causality conditions (nondegeneracy of local light cones) and a proper-time 
time foliation.^ For calculational purposes, these lattice configurations are then ro- 
tated to Euclidean signature and the path integral over this class can in principle 
be performed. Of course, since the physics one hopes to describe ultimately by 
these theories has Lorentzian character, one will have to perform an "inverse Wick 
rotation" back to Lorentzian spacetime eventually, never mind whether the compu- 
tation at an intermediate step took place in a purely Euclidean or in a Euclideanized 
Lorentzian framework. 

The simplest implementation of Euclidean DT based on the lattice Regge ver- 
sion of the Einstein-Hilbert action (the inclusion of a cosmological term being un- 
derstood) does not seem to lead to a theory with an interesting continuum limit. 
Even if this is the case, it is in principle possible that by adding more terms to 
the bare lattice action and suitably tuning the associated new coupling constants, 
an interesting continuum theory may emerge after all. This possibility has been 
investigated in the past [19] , as well as more recently [20] , but there is no conclusive 
evidence at this point that these modified Euclidean models can reproduce the phys- 
ical properties of quantum gravity from CDT, the Lorentzian lattice gravity theory 
to which we will turn next (see also [21] for a variety of reviews of the subject). 

Fig. 4 illustrates the general construction of a four- dimensional CDT triangu- 
lation. We take space to be compact and with the simplest topology, that of the 
three-sphere S^. In addition, we assume a discrete proper-time foliation and rep- 
resent the spatial geometry at each integer proper time t by a three-dimensional 
simplicial manifold, given as some configuration of Euclidean DT in terms of equi- 
lateral tetrahedra. By assumption, the tetrahcdra are flat in the interior, which 
means that their geometric properties are uniquely specified by the length of their 
edges, which is some number > (the same for all edges). To obtain a four- 
dimensional Lorentzian simplicial manifold with signature ( — h ++), we still must 
fill in all intervals [t, t + 1] between consecutive spatial slices. This can be done by 
using two types of geometrically distinct four-simplices, which again by assumption 
are flat in the interior, but this time with Lorentzian signature. The two different 
types are the (4,1)- and the (3,2)-simplex depicted in Fig. 4, together with their 
time-reversed counterparts. The (4,l)-simplex has as its "base" one of the spatial 
tetrahedra contained in the triangulated constant-time slice. (The "4" in the la- 
bel (4,1) refers to the four vertices contained in slice t that span this tetrahedron; 
similarly, the "1" refers to the single vertex shared with slice t + 1. An analogous 
labehng has been used for the (3,2)-simplex.) All that remains to be done to flx 
the geometry of the four-simplices is to assign lengths to the edges that have their 
end points in adjacent slices, and whose time labels therefore differ by one unit. We 
choose them to be all time-like and of equal (absolute) length at > 0, which in our 

®Note that there is no strict physical requirement that individual path integral histories must 
be causal; individual histories are not physical, observable quantities, only expectation values 
computed in the ensemble of histories are. 
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Figure 4: A triangulation in CDT consists of four- dimensional triangulated layers 
assembled from (4,1)- and (3,2)-simplices, interpolating between adjacent integer 
constant-time slices (left), which in turn are triangulations of 5*^ in terms of equi- 
lateral tetrahedra. Each purely spatial tetrahedron at time t forms the interface 
between two (4,l)-simplices, one in the interval [t — 1, 1], and the other in [t,t + 1], 
as illustrated on the right. Although a (3,2)-simplex shares none of the five tetra- 
hedra on its surface with a constant-time slice (the tetrahedra are all Lorentzian), 
it is nevertheless needed in addition to the (4,l)-building block to obtain simplicial 
manifolds with a well-defined causal structure. 

signature convention implies that their squared edge length is given by —afJ 

Our choice of causal geometries and length assignments has the added benefit 
that we can define a map that uniquely maps each Lorentzian CDT history to a 
Euclidean DT history. Let us start by parametrizing the relative length of the two 
lattice parameters and at by a positive real number a defined by a := —af/a^. 
Performing a rotation a — —a in the complex lower-half plane can be interpreted as 
changing all time-like length assignments of lattice links to space-like ones according 
to 

% = —aa^ — )■ = aa^. (18) 

In order that the Euclidean four-simplices obtained after this rotation satisfy triangle 
inequalities we require a > 7/12. The resulting triangulation represents a piecewise 
linear manifold with Euclidean signature. If one writes the Lorentzian Regge action 
as a function of a single lattice parameter a := and of a, the action behaves under 
the rotation (18) as one would expect naively from a rotation from Lorentzian to 
Euclidean spacetime, namely, 

iSL[a] = -SE[-a]. (19) 

The prescription (18) leading to (19) is the "Wick rotation" we had in mind in our 

''Note that at gives us an approximate distance measure between adjacent spatial slices labelled 
by integer-t, where the distance of a point in slice i + 1 to slice t is defined as the length of the 
longest geodesic from the point to the slice. 
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earlier discussion in Sec. 2.2. It transforms the original Lorentzian path integral 
with complex weights e^^^iT) with real weights e~'^^(^\ where by shght abuse 

of notation we use the same symbol T to denote the initial triangulation (with 
Lorentzian edge length assignments) and the one after rotation (which has identical 
connectivity, but purely Euclidean edge length assignments). Modulo the sign flip for 
the length assignments, the domain of the Euclideanized path integral is the same set 
T = {T} of triangulations as that of the original Lorentzian path integral. The set 
T is of course smaller than the set of all Euclidean triangulations one would obtain 
by gluing together the same Euclideanized building blocks, because it still carries 
an imprint of the causality conditions imposed on the Lorentzian triangulations. 

The fact that in DT and CDT we use standardized building blocks to construct 
the triangulations means that the Regge action takes on a very simple functional 
form. For the special case |a| = 1 we have after the Wick rotation only a single type 
of building block, the equilateral four-simplex with all link lengths equal to a = a^. 
The Regge form of the Einstein- Hilbert action becomes 

SE[-a = -1; T] = -noNo{T) + n^N^iT), (20) 

as is well-known from Euclidean DT quantum gravity. In (20), Nq(T) denotes 
the number of vertices in the triangulation T, and Ni{T) the number of its four- 
simplices. The coupling hq is related to the gravitational coupling constant G via 
I/kq oc Ga^, and should be identified with a^K/G, where A is the cosmological 
constant. 

Whenever |a| ^ 1, we retain the two different building blocks (of type (4,1) 
and (3,2)) after the rotation, and the action will depend on their total numbers, 
N^f'^^ and Nf''^\ separately instead of only on their sum A^4 = Nf'^'^-\-Nf''^\ It is 
convenient to parametrize the resulting Euclideanized Regge action in the form 

SE[-a-T] = -{Ko + QA)No{T) + K,(^Ni'''\T) + Ni'''\T)^+ (21) 

A(^Arf'2)(r) + 2Arf'i)(r)), 

where the asymmetry parameter A is a function of a such that A(a = l) = 0. 

We note that A appears in (21) on a par with the other two coupling constants, 
Kq and K4. In what follows, we will treat it as a third independent coupling constant. 
The reason for doing this - despite the fact that it has no immediate interpretation in 
the Einstein-Hilbert action - is that in the region of phase space (the space spanned 
by the three couplings Kq, and A) where we observe interesting, apparently 
continuum physics, the entropy of geometries is as important as the contributions 
coming from the bare action term. To make this more explicit, one can rewrite 
the Euclidean partition function of the theory as a sum over the counting variables 
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Np'', Np'> and A^o according to 

Z{ko,k,,A) = ^e-^-[^l (22) 

T 

where J\f{N[^'^\ N^^''^\Nq) is the number of triangulations with N^^'^^ four-simphces 
of type (4,1), AQ ' four-simphces of type (3,2) and A^o vertices. Introducing the 
notation ci = Nq/N^'^'^^ and C2 = A^4^'^^/iV4^'^'', the leading-order behaviour of this 
combinatorial quantity in the large-volume limit is known to be of the form 

J\f{Ni'''\ iVf iVo) = e^(^i''^^)<'''+-'-, (23) 

where "s.l." denotes subleading terms in n['^'^\ and the q typically have some 
boundedness properties. Since in the same limit the action (21) can be similarly 
approximated by Se = fici,C2)N[^'^^ + s.l., it implies that in the region of phase 
space where the four-volume can become large, both J\f and e"*^^ have the same 
functional form and are potentially of the same magnitude. It turns out that this 
is the same region where we observe interesting continuum-like physics. Because of 
contributions from both "energy" and "entropy" , it is clear therefore that the effec- 
tive action governing physics in this non-perturbative region can be very different 
from the "naive" Einstein-Hilbert action, justifying our inclusion of A as a tunable 
parameter in the bare action. 

To summarize: taking as our starting point spacetimes with Lorentzian signa- 
ture, we can consider the transition amplitude between an initial and a final spatial 
three-geometry, [gl ] and [g^ ] separated by a proper time t. We can then regular- 
ize the theory, using CDT, representing three- geometries by equilateral Euclidean 
triangulations and spacetime geometries by causal, Lorentzian triangulations with 
a discrete proper-time foliation. In the CDT framework, each of the latter can be 
rotated to Euclidean signature, leading to a regularized, Euclideanized sum-over- 
histories. What remains to be done is to "remove the regulator", that is, take the 
lattice spacing a to zero. Denoting the initial and final spatial triangulations by 
and T^^\ we thus arrive at the prescription 

GE([9i% [9i%t, Ko, K,, A) lim J] e-^-1^1, (24) 

which can be viewed as the four-dimensional generalization of the two-dimensional 
loop-loop amplitude G{£i,£2,t) introduced in (15) above. For a more detailed de- 
scription of the CDT construction we refer the interested reader to [22, 23, 24]. 
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3 The phase diagram 



Contrary to the situation in two dimensions, we cannot calculate the amplitude 
(24) analytically. However, we can extract a lot of non-trivial, non-perturbative 
information by performing Monte Carlo computer simulations . This will usually 
start with an investigation of the structure of the space of coupling constants (the 
"phase space" of the underlying statistical system), in particular, trying to identify 
regions associated with a second-order phase transition, where according to standard 
lore one can hope to obtain continuum physics . 

Let us highlight two technical aspects related to our implementation of the com- 
puter simulations. Firstly, rather than fixing specific boundary three-geometries T^^^ 
at times and t, we take time to be periodic. Although this is strictly speaking 
in contradiction with imposing causality (it introduces closed time-like curves), in 
practice it turns out to not affect results. The nature of the ground states of ge- 
ometry is such that by choosing t sufficiently large - assumed from now on - the 
boundary condition becomes irrelevant. 

Secondly, as we have discussed, the action (21) depends on three coupling con- 
stants, one of which, k,4, can be identified with the cosmological coupling constant, 
multiplying the spacetime volume V in the action. In the computer simulations 
it is convenient to keep this four-volume fixed, which means that the cosmological 
constant does not really play a role. We compensate for this by performing separate 
simulations at different (fixed) spacetime volumes. From these wc can in princi- 
ple reconstruct results which depend on the cosmological constant via a Laplace 
transformation. 



We are therefore left with two coupling constants, Kq and A. The corresponding 
phase diagram is shown in Fig. 5 [25] and exhibits three distinct phases, labelled A, 
B and C. Phase C appears to be the one relevant for contimmm physics, because only 
there do we observe extended four- dimensional universes [26]. A careful numerical 
analysis reveals strong evidence that the transition between phases C and A is first 
order, whereas between phases C and B we find a second-order transition [27] . This 
very exciting result implies that the B-C phase transition line is a candidate for a 
region in the coupling-constant plane where genuine UV continuum limits may exist, 
defined by approaching specific points on the line. Conversely, moving away from 
the transition line into phase C corresponds to going towards an IR limit. 



The reason why phase C is related to extended four- dimensional spacetimes is illus- 
trated in Fig. 6, which shows both a sample path-integral configuration generated 
by the computer during the Monte Carlo simulations, as well as the associated quan- 
tum observable, obtained by averaging in the ensemble. While of course we have 




(25) 
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Figure 5: The phase diagram of CDT quantum gravity in the {kq, A)-plane. 



access to the complete geometric information of the quantum spacetimes that are 
generated, only a single degree of freedom is depicted here, the three-volume of a 
spatial slice of the quantum spacetime as a function of proper time. The time ex- 
tension in a given simulation is always fixed (in the case at hand to 80 discrete time 
steps). What we observe in Fig. 6 is that the universe does not make use of the full 
time interval available, but has a non- vanishing volume only on a connected subset 
of the time axis.^ 

A quantitative piece of evidence in favour of a /onr-dimensional extended uni- 
verse is the fact that its time extension (not counting the stalk) scales like nI^'^ 
when the total discrete four-volume A^4 of the universe used in the simulations is 
varied. Similarly, its discrete three- volume A^3(t) scales like N^^'^. Contrary to one's 
naive expectations, these findings are highly non-trivial, because they have been 
derived in a non-perturbative, background-independent path integral formulation. 
The simplicial building blocks of our regularization are four-dimensional, but since 
assembling them is only dictated by the Boltzmann weight e~^'^^'^^ without any refer- 
ence to a four- dimensional background, there is no reason why the resulting object, 
extrapolated to infinite lattice volume, should be four- dimensional on any scale. 

This is specifically true in the non-perturbative regions of phase space where the 
entropic contributions to the effective action compete with those coming from the 
classical bare action, as explained above. In these regions it can easily happen that 
a type of configuration is entropically favoured that has no resemblance at all with 

^Since we impose the kinematical constraint that the spatial volume at fixed t cannot become 
smaller than 5 tetrahedra - the minimal number required to build a simplicial manifold of topology 

- the volume never vanishes completely. More precisely, what we observe in addition to the 
bell-shaped part of the volume profile is the formation of a distinct "stalk" which is close to the 
minimal size of 5 everywhere. 
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Figure 6: The three- volume of spatial slices as a function of proper time in phase C. 
Shown are a sample configuration of the volume profile, as well as the expectation 
value of the same quantity. 

an extended four-dimensional universe. Just from looking at the volume profiles, it 
is obvious that something like this does indeed happen in phases A and B, which 
as a result do not appear to have any classical limit resembling general relativity 
[26]. However, even in phase C the observed quantum universe is truly an outcome 
of non-perturbative dynamics, not a consequence of the dominance of the classical 
action.^ 

The fact that the path-integral measure can play a crucial role in determining the 
non-perturbative dynamics was a main lesson learned already earlier in the context of 
four- dimensional DT quantum gravity. When one considers a path integral ensemble 
of geometries obtained from gluing four-dimensional equilateral Euclidean simplices, 
with the only constraint that the topology should be that of S*^, one ends up with 
a universe of vanishing linear extension and infinite Hausdorff dimension [30]. This 
makes the situation depicted in Fig. 6 all the more remarkable! 

3.2 The effective action 

However, the surprises do not stop here. The smooth curve in Fig. 6 represents 
the expectation value of the volume profile, that is, the average over path integral 
configurations measured in the Monte Carlo simulations. For A^4 sufficiently large 
this curve is very precisely fitted by the function 



^Since we are working in Euclidean signature, dominance of the classical action would be fatal 
for the path integral, because of the action's unboundedness from below. In phase C, this instability 
is cured by the entropy of "microstates" or, in other words, the path-integral measure [28, 29]. 




(26) 
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where i denotes (integer) lattice time, the total number of four-simplices and 
N3{i) the number of tetrahedra at time i [31, 32], and Sq is a constant. 

Can the functional form of the expectation value found in (26) be obtained di- 
rectly from an action principle? The answer is yes [28] . A long time ago, Hartle and 
Hawking explored a minisuperspace approach to quantum gravity , where all gravi- 
tational (field) degrees of freedom at a fixed time are represented by a single number, 
the so-called scale factor or, equivalently, the total three- volume of the universe. 
Taking this classically reduced formulation as the starting point of the quantization, 
finding a quantum theory of gravity is reduced to a quantum mechanical problem 
in one variable, the scale factor a{t) [33]. 

The volume profile (26) of the emergent extended universe found in phase C 
of CDT quantum gravity can be derived from an "effective" action for the three- 
volume, namely, 

where t denotes proper time, k2 is a numerical constant and A is a Lagrange mul- 
tiplier, not a cosmological constant, because the total four-volume V4 is kept fixed 
in the simulations. Intriguingly, one obtains exactly the same expression (up to 
an overall sign) when plugging a spatially homogeneous and isotropic ansatz for 
the metric g^u{x) into the Euclidean Einstein- Hilbert action, and re-expressing the 
dependence on the scale factor in terms of the three- volume V^{t) oc a?{t). The 
solution to the equations of motion derived from (27) is the Euclidean de Sitter 
universe (a round four-sphere), which as a function of proper time t results in the 
cos^(t/Ki^''^)-dcpcndence of eq. (26). 

Despite the fact that they lead to very similar results for the dynamics of the 
scale factor, let us stress that conceptually there is a big difference between the 
ansatz of Hartle and Hawking, who simply assumed a minisuperspace reduction 
from the outset, and studying the effective dynamics of (the expectation value of) 
the scale factor in a full theory of quantum gravity, as we are doing. The only 
small but important reminder of the non-perturbative origin of the action (27) is 
its overall sign, which is opposite to that found in Euclidean cosmology. It can be 
attributed directly to "entropic" contributions to the effective action. The solutions 
to the equations of motion are of course not affected by this sign difference. A 
discretization of the effective action (27) has the functional form 

Sdiscr = h2_^ y j;^^ h k2N^' [t) - XNslt) j . (28) 

We have managed to reconstruct it in detail from the simulation data for the volume- 
volume correlator (V3(i)V3(f )), and have also shown that the quantum fiuctuations 

-•^^Thc formula is of course not valid in the stalk, where N^ii) ~ 5. 

^^This rather crude approximation is borrowed from standard cosmology, where homogeneity 
and isotropy are assumed to give a realistic description of our universe on the very largest scales. 
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around the de Sitter "background geometry" are well described by the action (28), 
yet another non-trivial result [32] . 

The same data have allowed us to relate the continuum coupling constant G in 
(27) to the constant ki in (28) according to 



where a is the lattice spacing and C4 is essentially the volume of a four-simplex (for 
lattice spacing a = 1), but depends weakly on the ratio between n'^'^^ and A^,^'''^'' 
(since the (4,1)- and (3,2)-simplices only have identical four-volumes when a = 1). 
This ratio, as well as the value of the constant Sq, defined in eq. (26), depend on the 
choice of the bare coupling constants Kq and A in phase C. 

Let us consider a typical choice for these couplings, (ko,A) = (2.2, 0.6), position- 
ing us in the interior of phase C. At this point in phase space, we have measured ki 
and with the help of (29) expressed Newton's constant and the Planck length £p in 
terms of the lattice spacing, resulting in 



From the identification of spacetime with a Euclidean de Sitter universe we have that 
V4 = Stt'^R^/S = C/^N/^a^, where C4 is the same quantity that appeared in (29). For 
the range of four- volumes used in the simulations, A'4 e [45.000, 360.000], the linear 
size ttR of the quantum de Sitter universes lies between 12 and 21 Planck lengths 
Ip. The small size of our universes is compatible with the fact that the observed 
quantum fluctuations in the three-volume are quite substantial, as illustrated by 
Fig. 6 (see also Fig. 7). For larger universes, the volume fluctuations will quickly 
become irrelevant. 

However, in order to investigate quantum properties of spacetime at Planckian 
and even sub-Planckian length scales, we want to do the opposite, namely, make the 
universes smaller and in this way increase the small-scale resolution of the simula- 
tions. How can we improve on (30) such that a single Planck length ip corresponds 
not to just half a lattice spacing, but to many lattice spacings a? From eqs. (29) and 
(30) it is clear that when ki goes to zero, £p can become much larger than a. The 
question is whether we can adjust k\ to go to zero. Since ki depends on the bare 
coupling constants kq and A, we have performed a scan of phase C to determine its 
qualitative behaviour [32]. Moving toward the A-C phase transition, ki is indeed 
decreasing, without going all the way to zero in the range of coupling constants 
scanned so far. Approaching the B-C phase transition is more difficult, because the 
system undergoes a second-order transition, and we observe a corresponding critical 
slowing-down. As far as we can tell from the numerical data at this stage, ki does 
not decrease when we approach this transition. However, as we will see in the next 
section, having ki go to a fixed value different from zero is actually the behaviour 
predicted at an ultraviolet second-order transition line, and therefore compatible 
with the continuum scenario we have appealed to earlier. 




(29) 



ki 3^6 ' 





(30) 
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3.3 Making contact with asymptotic safety 

Let us return to the renormalization group equation (2), which was formulated in 
terms of the dimensionless couphng constant G — G E'^. Now that we have a UV 
cut-off, the lattice link length a, we can instead form the dimensionless quantity 
G — Gjo? . From (29) it can essentially be identified with the inverse of fci, which 
we can measure. Wc can reformulate the renormalization group in terms of the new 
short-distance cut-off as 

An 

G[a) = c?G{a), a— ^ -13(G), /3{G) = 2G - cG'' + ■ ■ ■ , (31) 

where c depends on the constant u of eq. (2). Near the putative non-Gaussian UV 
fixed point G*, we can expand G and ki to lowest order in a according to 

G(a) = G* - Ka^ h(a) ^ kl + Ka\ (32) 

for some K, K, where the approach to the fixed point is governed by the exponent 

c^-P'{G*). (33) 

As explained in Sec. 2.1, in standard lattice theory one would now relate the lattice 
spacing near the fixed point to the bare coupling constants with the help of some 
correlation length ^. However, in four- dimensional quantum gravity we do not yet 
have a suitable correlation length at our disposal which could play this role. 

In search of an alternative, let us first consider the equation V4 = N^a'^, which 
defines the dimensionful continuum four-volume V4 in terms of the number A'4 of 
four-simplices and the lattice spacing. If we could consider V4 as fixed, we could 
replace the a-dependence of (32) by a iV4- dependence, with the advantage that N4 
is a parameter we can straightforwardly control. Re-expressing eq. (32) in terms of 
A^4 yields 

k^{N,) = kl-K'N-'/\ (34) 

for some K'. Since we can measure ki, we could determine the flow to the fixed 
point. The question is now which lattice measurements we should perform in order 
to make eq. (34) applicable. Increasing A^4 while staying at a specific point (kq, A) 
in phase C does not correspond to keeping V4 fixed, because during this process 
the size of the quantum fluctuations in the three-volume decreases relative to the 
expectation value of the three-volume. (More precisely, we already know that the 
ratio goes to zero like l/N^^.) Conversely, if "physics" is to be constant, which 
includes a constant V4, that same ratio should also remain constant. 

We will use this observation as our definition for what we mean by a "path of 
constant physics". If wc had a correlation length available, we could increase A^4 
and simultaneously change the bare coupling constants in such a way that the ratio 
of the correlation length to the linear extension of the universe of volume N4 (both 
in terms of lattice units) stayed constant. In the absence of a suitable correlation 
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Figure 7: Left: Three- volume profile for given A^4, for specific values (kq, A) of 
the bare coupling constants. Also indicated is the magnitude of the three-volume 
fluctuations around the mean value. While the expectation value of the three- 

3/4 1/2 

volume scales like , the fluctuations only scale like . Right: identifying 
a path of "constant physics" in the kq-A plane. Starting at some point in phase 
C, a path moving toward the UV phase transition is created by increasing A^4 and 
simultaneously adjusting and A, such that the ratio of the size of the three- volume 
fluctuations and the expectation value of the three- volume remains constant. 

length, we will use the magnitude of the three-volume fluctuations instead, and 
identify a "path of constant physics" as a trajectory in phase C along which the 
discrete four-volume A^4 grows, but the accompanying change in the bare couplings 
Kq and A ensures that the three-volume fluctuations likewise increase, in such a 
way that the ratio between the magnitude of the fluctuations and the mean three- 
volume stays the same. Fixing this ratio forces us to change bare coupling constants 
when we increase iV4, in this way tracing out a path that moves toward one of the 
phase transitions bordering phase C, see Fig. 7 (right) for a schematic illustration. 
Preliminary results from computer simulations to determine the flow defined in this 
way indicate that it should start quite close to the B-C phase transition if it should 
resemble the flow line of constant physics shown in the figure, raising again the issue 
of critical slowing-down near the B-C line. 

4 Relation to Hofava-Lifshitz gravity 

As described above, our CDT data in phase C can be fitted well to the functional 
form (28), which in turn can be seen as a discretized version of the minisuperspace 
action (27). There is a residual ambiguity in the interpretation of the discrete time 
coordinate appearing in the identification (26), which can be thought of as an overall, 
finite scaling between the time and spatial directions. As we have emphasized, due to 
the entropic nature of the effective action, there is no compelling reason to take the 
geometric length assignments of the regularized theory literally. We have identified 
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the time "coordinate" t with continuum proper time in such a way that we obtain 
a round four-sphere, which is a perfectly legitimate and physically well-motivated 
choice. However, as we vary the bare couplings kq and A, the overall shape of the 
computer-generated universe changes in terms of the number of lattice spacings in 
the time direction relative to those in the spatial directions. Although this change 
is qualitatively in agreement with the change of a as a function of kq and A, there 
is no detailed quantitative agreement. 

Instead of choosing continuum time to be consistent with a continuum S^- 
geometry as one moves in phase space, one may be able to find a modified action 
which describes the observed behaviour without performing an overall time rescal- 
ing which depends on kq and A. This may be especially appropriate in the vicinity 
of the phase transition, where the length scales one is probing become increasingly 
Planckian, and one would expect significant contributions to the effective dynam- 
ics from terms not contained in the infrared form of the Einstein-Hilbert action 
including higher-order curvature terms. 

We will consider yet another generalization, which suggests itself because of the 
built-in anisotropy between time and space of the CDT set-up, namely, a defor- 
mation a la Hofava-Lifshitz [34]. A corresponding effective Euclidean continuum 
action, including measure contributions, and expressed in terms of standard metric 
variables could be of the form 

Sh = j ^^x dt N^[{Ki^K'^ - XK'') + (-7i?(3) + 2A + T/(^,,)) , (35) 

where Kij denotes the extrinsic curvature and gij the three-metric of the spatial 
slices, i?*^^-* the corresponding three-dimensional scalar curvature, the lapse func- 
tion, and finally V{gij) a "potential" which in Hofava's continuum formulation would 
contain higher orders of spatial derivatives, potentially rendering Sh renormalizable. 
In our case we are not committed to any particular choice of potential V{gij), since 
we are not imposing renormalizability of the theory in any conventional sense. 

An effective V{gij) could be generated by entropy, i.e. by the measure, and may 
not relate to any discussion of the theory being renormalizable. The kinetic term 
depending on the extrinsic curvature is the most general such term which is at 
most second order in time derivatives and consistent with spatial diffeomorphism 
invariance. The parameter A appears in the (generalized) DeWitt metric, which 
defines an ultralocal metric on the classical space of all three- metrics^^, and the 
parameter 7 can be related to a relative scaling between time and spatial directions. 
Setting A = 7 = 1 and V — Q'm (35) we recover the standard (Euchdean) Einstein- 
Hilbert action. 

^^The value of A governs the signature of the generahzed DeWitt metric 

which is positive definite for A < 1/3, indefinite for A = 1/3 and negative definite for A > 1/3. The 
role of A in three-dimensional CDT quantum gravity has been analyzed in detail in [35] . 
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Making a simple minisuperspace ansatz with compact spherical slices, which 
assumes homogeneity and isotropy of the spatial three- metric gij, and fixing the 

lapse to = 1, the Euclidean action (35) becomes a function of the scale factor a[t) 
(see also [36, 37, 38], as well as [39] for related work in 2+1 dimensions), that is, 

Srmm = j dt a{tf (^3(1 - 3A) ^ - 7 ^ + 2A + V{a)) . (36) 

The first three terms in the parentheses define the IR limit (which in Hof ava-Lifshitz 
gravity is assumed to include a fiowing of A to its "GR value"), while the potential 
term V{a) contains inverse powers of the scale factor a coming from possible higher- 
order spatial derivative terms. 

Our reconstruction of the effective action from the computer data is compatible 
with the functional form (36) of the minisuperspace action. If we were able to 
extract the constant k2 in front of the potential term in (28), it would enable us to 
fix the ratio (1 — 3A)/27 appearing in (36) [40]. At this stage, the precision of our 
measurements is insufficient to do so. The same is true for our attempts to determine 
V{a) for small values of the scale factor, which is important for understanding UV 
quantum corrections to the potential near a{t) — 0. Once we have developed a 
better computer algorithm which allows us to approach the B-C phase transition 
line more closely, investigating such Planckian properties and testing scenarios of 
Hofava-Lifshitz type will be within reach. 

4.1 Conclusions 

In constructing a theory of quantum gravity using Causal Dynamical Triangula- 
tions, one of our initial inputs was the Reggc action, which appears in the weights 
of individual spacetimes in the gravitational path integral. However, as we have 
emphasized repeatedly, the full effective action generated dynamically by perform- 
ing the non-perturbative sum over histories is only indirectly related to this "bare" 
action. Likewise, the coupling constant ki, which appears in front of the effective ac- 
tion and we view as related to the gravitational coupling constant G, has no obvious 
direct relation to the "bare" coupling kq appearing in the Regge action. 

Nevertheless, the leading terms in the effective action for the scale factor are 
precisely the ones present in (27) or, more generally, in the effective Horava-Lifshitz 
action (36), at least for sufficiently large values of the scale factor. The fact that 
a kinetic term quadratic in derivatives appears as the leading term in the effective 
action is perhaps less surprising, but that the correct powers of the (undifferenti- 
ated) variable N^li) appear in both the kinetic and potential terms in (28) is rather 
remarkable and very encouraging for the entire CDT quantization program. 

For the range of bare coupling constants and four- volumes investigated until now 
our results are compatible with the Einstcin-Hilbcrt action. Better data and more 
observables will be required to discriminate between a "pure gravity" behaviour and 
an anisotropic deformation a la Hofava-Lifshitz in the deep ultraviolet. A beautiful 
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feature of CDT quantum gravity is that entirely non-perturbative questions of this 
kind can be formulated explicitly and addressed with the non-perturbative lattice 
tools available, and - if one is lucky - be answered quantitatively. 
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